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Abstract 

We study constraints on the natural inflation model from the cosmic microwave 
background radiation (CMBR). Inflaton cj) for the natural inflation has a potential of 
the form V = A 4 [l — cos(</>/ \/2/</>)], which is parametrized by two parameters and 
A. Various cosmological quantities, like the primordial curvature perturbation and 
the CMBR anisotropy, are determined as functions of these two parameters. Using 
recent observations of the CMBR anisotropy by BOOMERANG and MAXIMA (as 
well as those from COBE), constraints on the parameters fs and A are derived. The 
model with lower than 8.5 x 10 18 (5.4 x 10 18 , 4.5 x 10 18 ) GeV predicts a power 
spectrum with index ns smaller than 0.95 (0.9, 0.85) which suppresses the CMBR 
anisotropy for smaller angular scale. With such a small ns, height of the second 
acoustic peak can become significantly lower than the case of the scale-invariant 
Harrison-Zeldovich spectrum. 



Inflation [fl] plays a very important role in modern cosmology. It provides a natural 
solution to the horizon and flatness problems. In addition, quantum fluctuation during the 
inflation becomes the origin of the density perturbation of the universe. Such a fluctuation 
also generates temperature perturbation in the cosmic microwave background radiation 
(CMBR) which is being measured very accurately by the satellite and balloon experiments. 

In the case of the "slow-roll inflation" |2|, [|, where the inflation is due to the energy 
density of a slowly-rolling scalar field, the scalar field, which is called "inflaton," has to 
have a very flat potential. It is, however, difficult to guarantee the flatness of the potential 
because radiative corrections to the scalar mass are in general quadratically divergent. 
Thus the most natural value of the scalar mass is as large as the cutoff scale which is 
naturally the Planck scale. With such a large mass parameter, it is impossible to build a 
realistic model of inflation. 

One way to obtain a very flat potential is to consider the "natural inflation" 0, [IJ where 
a (pseudo-) Nambu-Goldstone (NG) field is used as the inflaton.!*]] If a global Abelian 



symmetry, which we call U(l)x, is spontaneously broken, massless scalar field shows up 
due to the Nambu-Goldstone's theorem. If the U(l)x symmetry is softly broken, such a 
scalar field acquires a non- vanishing potential, and the height of the potential is controlled 
by soft-breaking parameters. In this case, hierarchy between the scalar mass and the 
Planck (or the cutoff) scale is naturally stabilized. 

In the natural inflation, the inflaton has a particular form of the potential: V = A 4 [l — 
cos(0/\/2/</))]- Therefore, observable quantities are quite predictive and it is interesting 
to study its consequences. Indeed, the density fluctuation from the natural inflation was 
already extensively discussed in connections with the large scale structure and the COBE 
observations 0. Recently, however, BOOMERANG and MAXIMA [0] reported more 
accurate observations of the CMBR anisotropy up to the multipole / ~ 0(100). With 
those new observations, we can have a better constraint on the natural inflation scenario. 
(For discussions on other inflation models, see Refs. || ||.) 

In this letter, we consider constraints on the natural inflation model from the CMBR. 
For an accurate estimation of the power spectrum, we numerically follow the evolution 
of the inflaton field. The power spectrum deviates from the scale-invariant Harrison- 
Zeldovich spectrum when is relatively small, and the CMBR anisotropy from the natural 
inflation can be significantly different from that from the scale-invariant power spectrum. 
Furthermore, taking the constraint on the height of the second acoustic peak seriously, we 
may obtain an upper bound on the scale fy. 

We consider a model where NG boson is used as the inflaton. As an example, let us 
start with a model with the following scalar potential 



V{X) = A (|X| 2 - v 2 ) + AMt N (X N + X Nj ) , (1) 

where A is a complex scalar field, and A is a positive integer. We assume that smallness 

^Another way to guarantee the flatness is to introduce supersymmetry. In the supersymmetric models, 
quadratic divergences are cancelled out among bosonic and fermionic loops. In this paper, we do not 
consider such a possibility. 
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of model parameters should be protected by some symmetry; thus, we take A ~ 0(1) 
and v of the order of the reduced Planck scale M* ~ 2.4 x 10 18 GeV. On the contrary, 
when A = 0, there is an Abelian symmetry to rotate the phase of X, which we call U(l)x 
symmetry. Smallness of the A parameter is protected by U(l)x- In the following, we 
consider the case with A <C 1. 

The U(l)x symmetry is broken by the vacuum expectation value of X. In this case, 
it is convenient to parameterize the X field as 

where a and are real scalar fields. The real component o acquires a mass as large as v 
and is irrelevant for our discussion. On the contrary, the imaginary component is the 
(pseudo-) NG mode and is massless if A = 0. When A ^ 0, acquires a potential asp^ 2 



V{4) = A 4 



1 — cos 




where we added a constant to the potential for the vanishing cosmological constant, and 

4 _ 2Av N v_ 
~ M N ~ A ' ~ N' 

As one can see, the height of the potential is controlled by the soft breaking parameter A 
and the potential can be very flat with a small value of A. Since the flatness of the potential 
is protected by the (approximate) symmetry, quantum correction does not destabilize the 
flatness. The potential given in Eq. @ is our starting point. 

V(4>) has a minimum at = (mod 2 \f2ix f^). Expanding the potential around the 
minimum, we obtain V = |m^,0 2 + O(0 4 ), where the mass of </> is given by 

A 2 

At this level, is stable. If <p couples to other fields, however, it may decay. For example, 
X may couple to a fermion ipQ with standard-model gauge quantum numbers with the 
interaction £ int = uqXi/jqiPq + h.c. (Here, yq is a coupling constant.) Once X acquires the 
vacuum expectation value, i/jq becomes massive and the process — > 4>q?Pq is kinemat- 
ically blocked. At the one-loop level, however, decays into the standard-model gauge 
boson pairs. The decay rate for this process is given by 

^ dim(Gi) (b%ai\ 2 ml 



E 



1 32vr \ 4vriV ) fi 



(6) 



# 2 There may be other terms with higher periodicity in general. If the potential of cf> is due to U(l)x 
breaking spurion, however, the potential of (j> is expected to be of the form of Eq. (||). For example, if A 
originates to a spurion with charge N, sub- leading term is ~ A 2 cos(v / 2</ ) //0 + a), where a is an unknown 
phase. For realistic natural inflation, 4<1 and the sub-leading contribution is negligible. 
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where oti is the coupling constant for the gauge group Gi = SU(3)c, 577(2)l, and U(1)y 



while bj\ the /3-function coefficient of ipg + tpg. Here, the sum is over the standard model 



gauge groups. In our calculation, we approximate the formula as = 10~ 5 m^/ /J. 

We use the potential (0) and identify the field as the inflaton. We started with 
a specific example of the potential ([I]). However, notice that the following discussion is 
independent of the structure of the underlying model; once the potential (|f) is given, the 
following results all hold. 

If the field is displaced from the minimum in the early universe, it rolls towards the 
minimum of the potential. When t^z T^ 1 , the scalar field obeys the following equation 
of motion: 

4> + 3H<j> + V' = 0, (7) 

where the "dot" is the derivative with respect to time t while V = dV/ dip. Here, H is the 
expansion rate H = a/a with a being the scale factor. 

When the displacement of from the minimum is larger than ~ M#, slow-roll condi- 
tions, e < 1 and are satisfied, where 

1 ( M*V'\ 2 MlV" 
V 



2 V V J 1 ' V 

Then, the energy density of the universe is dominated by the potential energy of the 
inflaton field, and the universe is in the de Sitter phase. In this case, d > and the 
comoving scale grows faster than the horizon. 

When becomes less than the Planck scale, the slow-roll conditions do not hold. Then, 
d becomes negative and the inflation ends. This happens when 

2 = V. (9) 

After this epoch, the scalar field starts to oscillate and the amplitude of the oscillation 
decreases. Eventually, the amplitude becomes so small that the scalar potential is well- 
approximated as V ~ |m^0 2 . Then the energy density of is proportional to a~ 3 (as far 
as the decay process is neglected). In this period, we can solve the Boltzmann equation 
for the energy density of instead of following the motion of 0, since the change of the 
energy density during the time scale m^ 1 becomes very small: 

+ ZHp^ = -r^. (10) 

In addition, the radiation energy density p ra( j obeys 

Prad + 4-ffp ra d = (11) 



When t ~ T^ 1 , the inflaton decays. Then the energy density of is converted to that of 
the radiation and the universe is reheated. The reheating temperature is approximately 
given by T R ~ JyIm*. 
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We follow the evolution of the inflaton field numerically. Before showing the numerical 
results, however, it is instructive to discuss the qualitative behavior of the result with 
some approximation. To make the situation clear, we consider the case where the initial 
value of the inflaton to be < < \[2ixf t ^ n although the results do not depend on this 
assumption. Solving Eq. (Jfty with the slow-roll approximation, we obtain the amplitude 
at the end of the inflation as 



-'end 



12) 



When > en d, a > and the universe inflates. The e-folding number during the inflation 
is estimated as 



In 



cos 



^cnd 



/2V2U 



cos 2 (0/2V / 2/<^ 



(13) 



The scale of the COBE observation (& C obe ^ 7.5a H JL0| with H = lOO/i km/sec/Mpc 
and ao being the present expansion rate and the scale factor, respectively) exits the horizon 
when N e ~ 50 — 60. From Eqs. ( |12"1) and (|13|), it is clear that, when 3> M*, our horizon 
scale is affected only by the inflaton dynamics with <C f^. In this case, in particular, the 
inflaton amplitude for the COBE scale 0cobe becomes much smaller than fy. Then, we 
can approximate the inflaton potential as V ~ |m 2 2 , and all the predictions are the same 
as those from the chaotic inflation model with a parabolic potential with mass parameter 
m«£. If fcf, ~ 0(M*), on the contrary, 0cobe ~ 0(f<j>) and V(<p) cannot be approximated 
as above. In this case, we may observe some peculiar signal due to the natural inflation 
potential (H). 

Quantum fluctuation of <fi during the inflation becomes the origin of the density fluctu- 
ation. At the time of the horizon exit, the power spectrum of the curvature perturbation 
is given by 



V 1/2 (k) 



H 2 



2vr|0| 



(14) 



k=aH 



V(k) is often approximated using the power index ris; denoting the curvature perturbation 
around k ~ k as V(k) ~ V(k) x [k/k) ns<yk \ the index is given by |T0[] 



M 2 

n s (k) = l-6e(k) + 2 V (k) ^1-^ 



l + cos 2 (0/2v^2/ ) 
1 - cos 2 (0/2 v^/*) 



(15) 



k=aH 



Thus, the spectrum deviates from the scale- invariant one if 0cobe ~ 0(f ( f > ), which is 
realized when ~ 0(M*). Importantly, the spectrum index becomes smaller than 1, 
and the fluctuation for smaller scale is more suppressed compared to the case of the scale- 
invariant power spectrum. 

Now, we show the results of our numerical calculation. In our calculation, we first 
follow the evolution of the inflaton. From the period of the inflation to the time when 
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V(4>) is well approximated by the parabolic potential, Eq. (|7|) is solved. After that period, 
we follow Eqs. (|10D and (|TT|) until the inflaton decays. Once the radiation-dominated 
universe is realized, we use the simple scaling law (i.e., a 3 s = const., with s being the 



entropy density) to obtain the normalization of the comoving momentum k 

Following the evolution of the inflaton field during the inflation, we also calculate 
the curvature perturbation V as a function of the comoving momentum k as well as the 
spectrum index n^. We compared V(k) from our numerical calculation with that with the 
power-law approximation pi^p^){k) = (&/A;coBE) ns(fccOBE) ^(fccoBE), and found tliat the 
power-law approximation is in a good agreement with the numerical result. 

The temperature fluctuation observed by COBE sets a constraint on the primordial 
curvature perturbation V(k). In order to discuss the perturbation for the COBE scale, it is 
convenient to use the parameter Sh, density perturbation at the time of the reentry to the 
horizon. For the matter dominated era (with sizable contribution from the cosmological 
constant), Sh is given by |T2[ 

**(*) = l^V^k), (16) 
5 ilo 

where the above expression is for the flat universe, and 

J(SJo) = 5 no f± + ^_ii + eV 1 . (it) 

y\ oj 2 o ^ 7Q 14Q 14Q -r o j k j 

5h becomes the source of the temperature perturbation and is constrained so that the 
COBE observations [O] are reproduced 



IM^cobe)| = (1.94 ±0.15) x 1(T 5 

xQq 0.785-0.05 in n exp J_ .95(n s - 1) - 0.169(n s - 1) 

This constrains the vs. A plane. For large f^, as we mentioned, all the physical quantities 
are determined by the inflaton mass m^. In this case, the COBE scale exits the horizon 
when 4> ~ 15M*, and P 1/2 (A; C obe) ^ Im^jM*. Using Eqs. (§), QTB|) , and (Pp, the best-fit 
value of A for large enough is given by A ~ 8 x 10 15 GeV x {f ( f,/M, f ) 1 ~{oi fi = 0.4. 
When becomes small, higher order terms in the potential become more important and 
the best-fit value of A becomes smaller than the above simple expression. 

COBE also sets a bound on the spectrum index. The index for the COBE scale is 
constrained as JTD] 



™ S (£:cobe) = 1±0.2. (19) 

In Fig. 0, we plot ns(fccoBE) as a function of f^. Here, for each f^, we used the best- 
fit value of A. |^| Using the constraint (|i~9D , we can see that the natural inflation with 
frh ~ 3.8 x 10 18 GeV is inconsistent with the COBE observation. 



# 3 As indicated in Eq. (|l5|), ns(fccoBE) is insensitive to A for fixed value of We numerically checked 
the validity of this statement. 
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Figure 1: ns(fccoBE) as a function of fy. The best-fit value of A for the COBE-scale 
normalization is used. 



Now, we are at a point to discuss the CMBR anisotropy for the Z-th multipole C\. 



which is defined as 15 



(AT(x, 7) AT(x, f)> = — J2(2l + 1)QP^ • f 



(20) 



where AT(x, 7) is the temperature fluctuation of the CMBR pointing to the direction 7, 
and the average is over the position x. Theoretically, G\ is calculated once the transfer 
function Ti(k) is known: 

We used the CMBfast package ]16[ to calculate the transfer function and obtained C\ 
using the curvature perturbation V from the numerical calculation. The cosmological 
parameters used in our calculations are listed in Table [I]. In particular, suggested from 



the recent studies of the cosmological constant [18], we consider a model of flat universe 
with non-vanishing cosmological constant. 

In Fig. |2|, we plot C\ for several values of with the best-fit value of A for the COBE- 
scale normalization. With such a choice of A, Ci for smaller I is almost independent of f^. 
For larger I, however, Ci is more suppressed for smaller f^. This is because the curvature 
perturbation for smaller has a smaller index parameter n§. Notice that, with the best- 
fit value of A for the COBE-scale normalization, C\ becomes almost independent of for 
L>10 19 GeV. 
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h 


m 


w 


Q |T| 




Center value 


0.65 


0.019 


0.4 


Error (1-cr) 


0.08 


0.002 


0.1 



Table 1: Cosmo logical parameters used in our calculation. Here f2& and Qq are density 
parameters of baryon and total matter, respectively. We consider flat universe with cos- 
mo logical constant. 



Using the observations of the CMBR anisotropy by COBE, BOOMERANG and MAX- 
IMA, we can constrain the fundamental parameters and A. For this purpose, we cal- 
culate x 2 as P* 



xu,v=z l T~3ff' A E - <22) 

Here, C? hs is the observational data from COBE, BOOMERANG and MAXIMA for the 
Z-th multipole and 0"z )O bs is its error. The sum is over I available from the three experiments. 
The data used in our analysis are listed in Table |2|. In addition, C z th (/^, A) is the theoretical 
prediction to the Z-th multipole as a function of the fundamental parameters and A. 
Furthermore, we take account of the uncertainties in the cosmological parameters h, Qbh 2 , 
and Qq- We varied these cosmological parameters within the 1-cr error and calculated 
the variations in (^(/^A). We identified them as systematic errors and added them in 
quadrature to evaluate of th . 

In Fig. H|, we show the constraint on the vs. A plane. We shaded the region with 
X 2 < 43 (which corresponds to 95 % C.L. allowed region for the x 2- statistics with 29 
degrees of freedom). For a fixed fy, we obtain upper and lower bounds on A. With 
being fixed, the primordial curvature perturbation is an increasing function of A, and 
hence A is required to have a relevant value to explain the observed size of the multipoles. 
Notice that, for large enough fy, the upper and lower bounds on A are proportional to 
f X J 2 . This is because the observable quantities depend only on = A 2 /y/2f^ in this 
region. 

Now, let us comment on the case with 3> M*. In this case, the index ns becomes 
close to 1 and the theoretical predictions on C\ for I ^ 400 become larger than the observed 
ones. However, the asymptotic value of n s for 3> M* is 0.96, not exactly equal to 1. 
Thus, the constraints from the multipoles around the second peak is not as severe as those 
for models with the scale-invariant Harrison- Zeldovich spectrum. In addition, since we 

our calculation of y 2 , we do not take account of the correlations among data points. Using the 
correlation matrices given by BOOMERANG and MAXIMA, we derived a constraint on the vs. A 
plane, and checked that the changes of the upper and lower bounds on A for fixed are less than a few 
%. We also neglect the calibration uncertainties of the BOOMERANG and MAXIMA data sets (20 % 



for BOOMERANG and 8 % for MAXIMA, 1-cr in Ci |20|). We found that the inclusion of the calibration 



uncertainties may change the upper and lower bounds on A by a few %. 
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lett 


C? bs (/iK 2 ) 


07 iObs (/UK ) 


Reference 


3.1 


784 


473.5 


COBE 


4.1 


1156 


440.9 


COBE 


5.6 


630 


291 


COBE 


8.0 


864.4 


224.5 


COBE 


10.9 


767.3 


230.2 


COBE 


14.3 


681.2 


246.7 


COBE 


19.4 


1089 


326 


COBE 


51 


1140 


280 


BOOMERANG 


101 


3110 


490 


BOOMERANG 


151 


4160 


540 


BOOMERANG 


201 


4700 


540 


BOOMERANG 


251 


4300 


460 


BOOMERANG 


301 


2640 


310 


BOOMERANG 


351 


1550 


220 


BOOMERANG 


401 


1310 


220 


BOOMERANG 


451 


1360 


250 


BOOMERANG 


501 


1440 


290 


BOOMERANG 


001 


1 ( ou 


Oil) 




601 


1540 


430 


BOOMERANG 


77 


2000 


595 


MAXIMA 


147 


2960 


615 


MAXIMA 


223 


6070 


970 


MAXIMA 


300 


3720 


580 


MAXIMA 


374 


2270 


365 


MAXIMA 


447 


1530 


290 


MAXIMA 


522 


2340 


405 


MAXIMA 


597 


1530 


360 


MAXIMA 


671 


1830 


465 


MAXIMA 


742 


2180 


660 


MAXIMA 



Table 2: Cf hs and a hohs used in our analysis, which are from COBE BOOMERANG 
@| and MAXIMA fj experiments. 
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1 1 1 1 1 1 1 1 1 1 1 

100 200 300 400 500 600 700 800 900 1000 

/ 

Figure 2: CMBR anisotropy for the l-th multipole. The vertical axis is C\ = [l(l + l)/2ir]Ci. 
Here, we take / = 10 x 10 18 GeV (solid), 8x 10 18 GeV (dashed), 6x 10 18 GeV (dotted), and 
4 x 10 18 GeV (dot-dashed), and the best-fit values of A for the COBE-scale normalization 
are used. 



used all the data points to evaluate x 2 an d required x 2 < 43, the discrepancies for / ^ 400 
are not statistically significant enough to exclude the parameter region with 3> M*. 
Thus, in Fig. |3|, no upper bound on is obtained. However, it is interesting to consider 
the constraint from the second peak. As an example, following Ref. ||, we identified 
the highest data points for I < 400 and 400 < / < 600 as the heights of the first and 
second peaks, (7 z lst and Cf nd , respectively. Then, we combined the two data samples from 
BOOMERANG and MAXIMA to obtain 

Cf nd /Cl Bt = 0.38 ± 0.06. (23) 

Requiring Cf nA /C} st < 0.50, we obtain /^7x 10 18 GeV. In the future, more accurate 
measurements of the CMBR anisotropy will be able to check the validity of this upper 
bound. 

Another independent constraint is available from the cluster abundance. In Fig. |3|, 
we plotted the contours of the constant as, where as is the amplitude of mass density 
fluctuations on the scale of 8/i _1 Mpc. The observed value of as is given by a 8 = (0.56 ± 
0.06)f2o°' 47 P3I-P3 ^ s one can see ' ^ n case without reionization, as is consistent with the 



above-mentioned value in the parameter region preferred by the CMBR anisotropy. 

So far, we have not discussed the effect of the reionization after the recombination. 
Its effect is well parameterized by the following two parameters: the optical depth r and 

# 5 Constraint on erg is insensitive to the index parameter ns p2], and we neglect its dependence on n§. 
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< 10 - 
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< 10 : 
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2 3 4 5 7 10 20 
f^(10 18 GeV) 

Figure 3: Constraint on the parameters and A. The shaded regions are for \ 2 — 43 for 
(a) t = (no reionization), (b) r = 0.2, and (c) r = 0.4. The dotted lines are contours of 
constant as (0.4, 0.6, 0.8, 1.0, 1.2, and 1.4, from below). 
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the red shift z- lon at the time of the reionization [23|, 21 1. Due to the reionization, O/'s 
with I ^> z\l^ are suppressed by the factor e _2r while those with small I are unchanged. 
We calculated x 2 with the reionization with r = 0.2 and 0.4. For a fixed value of r, we 
took several values of z [on of O(10), and checked that x 2 is almost independent of z [on . 
The results are shown in Fig. || Since the reionization effect reduces Ci, larger value of 
the primordial perturbation is needed to obtain the CMBR anisotropy consistent with the 
observations. Thus, the preferred value of A becomes larger as r increases, as shown in Fig. 
|3|. In addition, when r = and ns = 1, theoretical prediction for C\ with I ^ 400 becomes 
larger than observations if we adopt the primordial curvature perturbation suggested by 
COBE. In this case, CMBR anisotropy prefers a index ns smaller than 1. However, if 
the reionozatoin effect is sufficient, it suppresses C\ with large I and ns < 1 makes the fit 
worse. Thus, for large r, parameter region with small f^, where ns becomes significantly 
smaller than 1, is excluded. In addition, consistency with the constraint from the cluster 
abundance becomes worse for larger r, as shown in Fig. || 

Finally, we briefly discuss possible improvement of the constraints with future obser- 
vations. With MAP p5] and PLANCK experiments, much better observations of the 
CMBR anisotropy will be obtained. It has been pointed out that MAP and PLANCK will 
determine the index ns with 0(1 %) accuracy [p7[] . In the natural inflation model, ns is 



sensitive to / . For example, 5.2 x 10 18 GeV < U< 5.7 x 10 18 GeV gives 0.89 < n s < 0.91. 
In addition, C\ itself will be determined much more accurately, and hence the theoretical 
prediction on the CMBR anisotropy will be more directly compared with the observation. 
Thus MAP and PLANCK will provide much better constraint on the natural inflation 
model than the present one. 
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Research from the Ministry of Education, Science, Sports, and Culture of Japan (No. 
12047201). 

References 

[1] A.H. Guth, Phys. Rev. D23 (1981) 347. 

[2] A.D. Linde, Phys. Lett. B108 (1982) 389. 

[3] A. Albrecht and P.J. Steinhardt, Phys. Rev. Lett. 48 (1982) 1220. 

[4] K. Freese, J.A. Frieman and A.V. Olinto, Phys. Rev. Lett. 65 (1990) 3233. 

[5] F.C. Adams, J.R. Bond, K. Freese, J.A. Frieman and A.V. Olinto, Phys. Rev. D47 
(1993) 426. 

[6] P. de Bernardis et al, Nature 404 (2000) 955. 



11 



[7] 

[8] 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 
[15] 
[16] 

[17] 
[18] 
[19] 
[20] 
[21] 
[22] 

[23] 
[24] 
[25] 
[26] 
[27] 



S. Hanany et al., |astro-ph / 0005 123 . 
W.H. Kinney, A. Melchiorri and A. Riotto, |astro-ph/ 0007375 
L. Covi and D. Lyth, |astro-ph/0008T65 



See, for example, D.H. Lyth and A. Riotto, Phys. Rept. 314 (1999) 1. 

See, for example, E.W. Kolb and M.S. Turner, "The Early Universe" (Addison- 
Wesley, 1990). 

S.M. Carroll, W. Press and E.L. Turner, Annu. Rev. Astron. Astrophys. 30 (1992) 
499. 

C.L. Bennett et al, Astrophys. J. 464 (1996) LI. 
E.F. Bunn and M. White, Astrophys. J. 480 (1997) 6. 
See, for example, W. Hu, Ph. D Thesis (|astro-ph / 9 508 1 26|) . 



U. Seljak and M. Zaldarriaga, "CMBFAST: A Microwave Anisotropy Code" 
(|http:/ /www. sns.ias.edu/~matiasz/CMBFAST/cmbfast. html ). 



J.R. Primack, b,stro-ph/0007187 . 
M.S. Turner, |astro-ph/99(J3D5I 



M. Tegmark and A.J.S. Hamilton, |astro-ph/9702019 



A.H. Jaffe et al, |astro-ph/0007333 . 

P.T.P. Viana and A.R. Liddle, |astro-ph/9902245| . 



S.D. White, G. Efstathiou and C.S. Frenk, Mon. Not. Roy. Astron. Soc. 262 (1993) 
1023. 

W. Hu and M. White, Astrophys. J. 479 (1997) 568. 
L.M. Griffiths, D. Barbosa and A.R. Liddle, |astro-ph/98f2T25" . 
MAP home page ( http: / /map. gsfc.nasa.gov 



PLANCK home page ( http:/ /astro.estec.esa.nl/SA-gcncral/Projects/Planck|) . 



J.R. Bond, G. Efstathiou and M. Tegmark, Mon. Not. Roy. Astron. Soc. 291 (1997) 
L33. 



12 



